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SUMMARY ' 

A quasi-three-dimensional finite-element analysis was used to analyze 
the edge— stress problem in four— ply, composite laminates. The seven laminates 
that were considered belong to the laminate family [6/(6-90)] g , where . 

0 < 0 <90. Systematic convergence studies were made to explore the existence 
of stress singularities near the free edge. The present analysis appears to 
confirm the existence of stress singularities at the intersection of the inter- 
face and the free edge. The power of the stress singularity was the same for 
all seven laminates considered. 

INTRODUCTION 

Many composite laminates develop high interlaminar stresses near free 
edges due to mechanical loads, thermal loads, or hygroscopic effects.. The 
interlaminar stresses can lead to delamination and eventual failure of the 
laminate. Several numerical analyses have computed stress distributions in a 
finite— width laminate subjected to uniform axial strain. This case is often 
referred to as the straight-edge or the edge— stress problem. In reference 1, a 
finite-difference scheme was used to analyze the edge-stress problem of a [±45] g 
laminate. In reference 2, a three-dimensional finite-element method was used to 
analyze the edge effects in [±45] g and [90/0] g laminates. In reference 3, a 

*Joint Institute for Advancement of Flight Sciences, George Washington 
University. 

This work was performed, in part, under NASA Grant NSG-1449. 



2 


quasi-three-dimensional finite-element analysis was used to analyze [0/ 90] ^ , 

[90/0] , [±45] , and several quasi-isotropic laminates. These analyses used 
s s 

coarse mesh models because of limitations on computing storages and times. 
However, all showed that the interlaminar stresses had very localized peaks 
near the free edge. If these localized peaks were caused by a singularity, 
as suggested in references 1 and 3, singularity analyses could be'developed 
and used to study delaminations. 

The present study investigated the existence of 9tress singularities at a 
free edge in several four-ply laminates that belong to the laminate family 
[0/(0-9O)] g , 0 < 0 < 90. To this end, quasi-three-dimensional finite-element 

analysis was used. The analysis, which is similar to that of reference 3, 
used eight-noded isoparametric elements. 

The interlaminar stresses are presented as distributions through the 
thickness at the free edge, and along the interface between the plies. The 
results are compared and discussed for seven different laminates. Stress 
distributions calculated from the analyses were used in a procedure to 
investigate the existence of stress singularities. 


SYMBOLS 



h 
r ' 

u,v,w 


strength of singularity, MPa 

half-width of the straight-edge laminate, m 

Young's modulus in i-direction, MPa 

shear modulus, MPa 

ply thickness, m 

nondimensional distance from the free edge, r 
displacement functions, m 


(b y)/h 
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u,v,w 

x,y,z 

a 

e 

o 

6 


v u 


{a} 


displacements in x-, y-, and z-directions, respectively, m 
Cartesian coordinates, m 
power of the singularity 

uniform axial strain imposed on the straight-edge laminate, 

e = 0.001 

o 

outer ply angle, deg 
Poisson’s ratios 

Cartesian stresses, O , O , O , 


xy- 


yz 


a , MPa 


zx 


Subscripts: 

i, j i,j = 1,2,3 

1,2,3 longitudinal, transverse, and thickness directions of a unidirec- 

tional ply 


ANALYSIS 

This section describes the laminates analyzed, the quasi— three— dimensional 
formulation, and the finite-element idealizations used in this study. This 
section also discusses the procedure used to investigate the existence of a 
stress singularity. 


Laminates 

The laminates considered in the analyses were [0/90] g , [15/-75] s , 

[30/-60] s , [±45] s , [60/-30] g , [75/-l5] g , and [90/0] As previously mentioned, 
these seven laminates belong to the family [6/(6“90)] g , 0 < 0< 90. The top 

ply angle, • 0, is shown in figure 1(a). 

Each ply was idealized as a homogeneous, elastic orthotropic material 

with properties (refs. 1, 2, and 3), 
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E n = 137.9 GPa (20 x 10 psi) 


E 22 " E 33 = 14 * 48 GPa (2,1 X 10 psi) 


G 12 " G 23 " G 3i “ 5.86 GPa (°* 85 x 10 P Bi > 


V 12 = V 13 = V 23 = °‘ 21 


The subscripts 1, 2, and 3 correspond to the longitudinal, transverse, and 
thickness directions, respectively, of the individual ply. 


Quasi-Three-Dimensional Analysis 

Figure 1(a) shows a long, symmetric laminate loaded in the x-direction. 
The laminate has a width of 2b and has four plies, each of thickness h. 

Away from the ends, the displacements in any x = constant plane were assumed 
to be 


u(x,y,z) = £ q x + U(y,z) 
v(x,y,z) = V (y , z) 
w(x,y,z) = W(y,z) 


( 1 ) 


where e is a uniform axial strain, and U,V,W are functions of the coor- 
o 

dinates y and z alone (ref. 1). Furthermore, because of symmetries 
inherent in the layups considered, U, V, and W satisfy the following 
requirements: 


U(y»z) = -U(-y,-z) > 

V(y,z) = -V(-y,z) > 

W(y,z) « -W(y,-z) j 


( 2 ) 
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and 

U(y,0) = V(y,0) = |^ W(0,z) = 0 

The displacement field of equations (1) should satisfy the equations of 
equilibrium and the stress-free conditions on the edges, y. = ±b, and top and 

bottom surfaces, z = ±2h. . 

Because exact solutions are not available for the complex three-dimensional 
problem, approximate methods were used in references 1 through 3. The present 
study was based on a finite-element analysis, like that in reference 3, except 
eight-noded quadrilateral isoparametric elements were used rather than three- 
noded triangular elements. 

Finite-Element Idealizations 

Because of the symmetries in the problem, only one quadrant of an 

x = constant plane was idealized. The quadrant (0 < y < b and 0 < z < 2h) 

is shown as the shaded region in figure 1(b). The displacement functions U 

and V were prescribed as zero on the y = 0 line, and the displacement 

function W was prescribed as zero on the z = 0 line. 

The shaded region in figure 1(b) was idealized by eight-noded isoparametric 

elements in the rectangular-mesh models shown in figure 2 and the polar-mesh 

model in figure 3. The rectangular-mesh model had a rectangular arrangement 

throughout. In contrast, the polar model had a polar mesh near the free edge 
* » 

and rectangular arrangement elsewhere. The rectangular mesh was used to com- 
pute the stress distributions, and the polar mesh was used to investigate the 
stress singularities. 

To study the convergence of the stresses near the free edge, three rec- 
tangular meshes were used. The medium mesh in figure 2(b) was obtained by 
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subdividing each element of the coarse mesh (fig. 2(a)) into four elements. 
Similarly, the fine mesh in figure 2(c) was obtained by subdividing each ele- 
ment of the medium mesh into four elements. The coarse mesh had 135 nodes and 
36 elements, the medium mesh had 485 nodes and 144 elements, and the fine mesh 
had 1833 nodes and 576 elements. 

The polar mesh had the polar arrangement in the region (19h < y < 20h, 

0 < z < 2h) near the laminate free edge. As in the rectangular-mesh case, 
three meshes — coarse, medium, and fine — were used. Only the polar fine mesh 
is shown in figure 3. The three polar meshes had the same number of nodes and 
elements as their corresponding rectangular counterparts. Near the free edge, 
however, the polar fine mesh was an order of magnitude finer than the rectangu- 
lar fine mesh. 


Analysis of Singularities 

The stress distribution along a radial line from a singularity can be 
expressed as (refs. 4 and 5) 

a = A^ -0 + 0(r" a+1 ) (3) 

where r is the distance from the singularity, and a are constants, 

and 0 (x -° +i ) represents terras of the order r -0 ^ and higher. For small 
distances r', the singular term dominates and equation (3) can be approxi- 
mated by ' 



or 

log o = log A 1 - a log r 


(4) 
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Hence, a log a vs. log r plot would be a line with a slope of -a and a 
a-intercept of A^. The a and are the power and the strength of the 

singularity, respectively. 

Equation (4) provided the basis for a precedure to investigate the 
existence of stress singularities at the laminate edge. This equation was 
fitted to the computed stresses. If the equation fitted these stresses well 
and if the slope of the fitted line was negative, a stress singularity with 
power ot was indicated. This log-linear procedure was verified in refer- 
ence 4 by applying it to several well-known singular— stress problems in two- 
dimensional elasticity. 

The log— linear procedure works best where the singular term strongly 
dominates the other terms. Therefore, the procedure should be applied only in 
the regions very close to the singularity. 

RESULTS AND DISCUSSION 

First, distributions are presented for interlaminar normal stress 

and the shear stress a for the various [0 / (0—90) ] , 0 < 0 < 90, laminates. 

xz s 

Then, the results are investigated for the existence of stress singularities 
at the free edge. For convenience, the intersection of the interface and the 
free edge (z = h; y = b) will be called the interface corner. Similarly, the 
intersection of the midplane and the free edge (z = 0; y = b) will be called 
the midplane corner. 

The uniform axial strain, £ q , was arbitrarily set equal to 0.001 through- 
out this study. 

Stress Distributions 

The stresses with the steepest gradients near the free edge were the 

normal stress a and the shear stress a . The distributions for these 
z xz 
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stresses, through the thickness and along the interface, are presented aiid com- 
pared in this section. 

Interlaminar normal stress, ° z *” 

• [0/90] laminate: Figure 4(a) shows the 0 distribution through the 

s z 

thickness along the free edge (y = b) for the three rectangular meshes (coarse, 

medium, and fine). In this figure and all subsequent figures, the fine mesh 

results are represented by a curve through the data; only the value at the 

interface is shown as a discrete value (diamond symbol) . The coarse and medium 

mesh results are shown by circular and square symbols, respectively. The solid 

symbols indicate the stresses in the top ply. As shown on the figure, the 

values of 0 for the three meshes agree closely except near the interface 
z 

(z = h) . • At the interface the three meshes produced noticeably dif ferent values 

In addition, 0 is discontinuous across the interface; this discontinuity is 
z 

slightly larger with progressive mesh refinement. Hence, the stress state shown 
in figure 4(a) suggests that a singularity exists at the interface corner. In 
contrast, at the midplane (z = o) the 0 z values agree for the three meshes 
and show no evidence of a singularity at the midplane corner. 

Figure 4(b) shows the average 0^ from both plies plotted against a nor- 
malized distance from the edge. The 0 z results from the three meshes are in 
excellent agreement for (b - y)/h > 0.08. However, at the free edge, y = b, 
the computed values are noticeably different for the three meshes, with 

the fine mesh producing the largest value. Therefore, results in both 
figures 4(a) and 4(b) suggest that a singularity exists at the interface corner. 

[±45] laminate: Figure 5(a) shows the 0 distribution through the 
s z 

thickness at the free edge (y = b) . As in the [0/90 ] g laminate, the values of 

a for the three meshes agree closely everywhere except near the interface 
z 

(z = h), where a steep gradient exists. Figure 5(b) shows the distribution of 
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0 along the interface. As expected, the results from the three meshes agree 

except very near the free edge. Further, the computed value of 0 z at the 

free edge was larger for progressive mesh refinement. These results for the 

£±45] laminate also suggest a 0^ singularity at the interface corner. 

Other laminates: Figures 6 and 7 show the 0, distributions through the 

z 

thickness and along the interface for [15/-75] g , [75/-15] g , [30/-60] g , and 

[60/-30] , respectively. In these figures only the fine mesh results are 
s 

shown. The distributions of o z near the free edge have steep gradients for 

all the four laminates, indicating a possible singularity. 

Comparisons: The results in figures 4 through 7 show that, except for 

[0/90] and [90/0] laminates, the interlaminar normal stress, 0 z> is compres- 
s s 

sive at the interface corner. In contrast, at the midplane corner only lami- 
nates [75/-15] and [90/0] develop compressive 0 . 
s s 

Interlaminar shear stress, ° xz * - 

[+45] laminate: Figure 8(a) presents the through-the-thickness distribu- 

tions of 0 at the free edge obtained using the three rectangular meshes for 
xz 

the [±45] laminate. Figure 8(b) shows the 0. xz distribution along the inter- 
s 

face. For all the meshes 0 closely agrees everywhere except near the 
interface in figure 8(a) and free edge in figure 8(b). The gradients are steep 
in this region, again indicating a possible singularity for the [±45] g laminate. 

Other laminates: For the other laminates considered the © xz distributions 

were similar to the [+45] results and, therefore, are not shown. Also, the 
average 0 values at the interface for a laminate and its corresponding com- 
plementary laminate (for example, [15/— 7 5] g and [75/-15] g ) differed only in the 
third significant digit. 
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Comparisons: For all laminates considered the 0 distributions had 

xz , 

the same sign along the interface. The [15/-75] and [75/-15] developed the 

largest 0 very near the free edge compared to other laminates. Near the 
X z 

free edge, the 0 stresses are larger (by at least a factor of 3) than the 
X z 

corresponding stresses in each laminate, except for the [0/90] ^ and [90/0] g 

laminates for which G is identically zero, 

xz 

For all laminates considered, G is identically zero at the midplane 

xz 

because of symmetry. Therefore, analysis of singularities is confined to the 
neighborhood of the interface corner. 

Singularities 

Both the G and G distributions suggest that singularities exist at 
z xz 

the interface corner. To identify these singularities, both G z and G^ 

stress data may be used in the log-linear procedure outlined in the analysis 

section. However, for the present edge-stress case (generalized plane strain), 

if one stress is singular at a point, then all stresses are and they all have 

the same singularity power, a. This conclusion follows from the observation 

that these stresses have similar functional forms, derived from the same order 

derivatives of the same stress functions, as shown in reference 6. Therefore, 

either o or G may be used in the log-linear procedure to determine the 
Z X z 

power of the' singularity. 

To decide which stress was better suited for the log-linear procedure, 

known singularities were examined. The case of a center-cracked plate showed 

that the log-linear procedure gave best results when the predominant stress was 

used. For the present study, the O xz stress was predominant. As previously 

mentioned, the 0 values near the edge were at least three times as large as 
xz 

the 0 values, except for the [0/90] and [90/0] laminates. Further, as 
z s s 
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pointed out earlier, the procedure works best when the singular term in 

equation (3) strongly dominates the other terms. The a distributions had 

xz 

the same sign all along the interface; in contrast, the a distributions 

z 

changed sign near the edge. This suggests that the singular term for O 

xz 

is dominant over a wider range of data than for a . For these reasons 

a was used in the log-linear procedure. 

The results for the rectangular and polar meshes were indistinguishable 

for distances greater than 0.04h from the free edge. Very near the free 

edge, however, the polar mesh yielded a greater concentration of data and, 

therefore, was used to analyze stress singularities. 

[±45] laminate.- Figure 9 shows the log-log plots of the absolute value 
s 

of 0-, along the interface for a [±45] laminate, plotted against r. The 
xz s 

stresses obtained from the three polar meshes are plotted in this figure. As 
expected, all three meshes show a linear portion with a negative slope, indi- 
cating a singularity. Before calculating the power of the singularity, the 
linear region for each set of data was estimated visually. (The two points 
closest to the edge were neglected because they represent nodes for the ele- 
ment at the singularity, and this element cannot be expected to yield accurate 
stresses.) Stresses from the coarse, medium, and fine meshes appear to be 
linear in the ranges of 0.0125 < r < 0.07, 0.0035 < r < 0.015, and 

0.0015 < r < 0.01, respectively. In these ranges a linear least-square- 
fitting procedure was used to determine the power (a) and the strength of 

the singularity in the equation 

log a = log A- -a log r 
xz i 

The following values were found for the three meshes. 
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POLAR MESH 

COARSE 

MEDIUM 

FINE 

Mesh index, i 

1 

2 

3 

Power, a 

0.256 

o! 182 

0.170 • 

r- 

Percent change in a, 

a. - a.,, 

— (100) 

a. 

' - 

41 

7 : 

Strength A^, MPa 

5.702 

7.509 

7.950 


As shown in the table, by refining the coarse mesh to generate the medium 
mesh, the a changed by 41 percent. For the fine mesh, the additional refine- 
ment produced only a 7-percent change, which is only about one-sixth the 
42-percent change. Further mesh refinement (dividing each fine mesh element 
ijito four elements) would probably cause an additional 1— percent change 
(one-sixth of 7 percent). This 1-percent change is neglibibly small. Similar 
arguments can be made for the strength of the singularity, A^. Therefore, 
the present fine mesh results were treated as converged values. 

To check the accuracy of the present procedure, the [±45] g singularity 
was computed using the stresses along nine different radial lines of the 

polar mesh.'* The calculated power of the singularity cl should be the same 
for all such lines. Figure 10 shows a and the corresponding ^ values 
plotted against <J>, the angular position of each radial line. These values of 
cl vary slightly, from 0.195 at 4* = 0 to an interface value of 0.170 at 
<J> = tt/ 2. At the interface, the 0^ z stresses were largest, as shown by the 
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largest value at <j> = tt/ 2 in figure 10. As mentioned earlier, the pres- 

ent procedure yields best results when the predominant stresses are used. 
Because the a ^ stress was largest at the interface (<j> = tt/ 2), the a value 
of 0.170 is believed to be accurate. Also, the absolute value of the correla- 
tion coefficient in the least-square procedure was closest to unity (0.9996) 
for <J> = ir/2. - " 

Other laminates.- The O stress was also used to identify the singu- 
xz 

larities for the [15/-75] , [75/— 153 > [30/— 60] , and [60/— 30] laminates. 

s s s s 

Figure 11 shows the log-log plots for these four laminates. For comparison 

the previous results for the [±45]^ laminate are also included in this figure. 

As pointed out earlier, for a laminate and its complementary laminate 

the a results differed only in the third significant digit. Therefore, 
xz 

the power of the singularity for a laminate and its complement is apparently 
identical. Also note that log-linear plots in figure 11 are parallel to each 
other. Therefore, the power of the singularities for [15/— 75] , [75/-15] , 
[30/-60] s , and [60/— 30] s are identical to that for the [±45] g case, that is, 
0.170. However, the three parallel lines in figure 11 are separated, indicating 
different singularity strengths, A^. 

Figure 12 shows a and A^ for various [0/ (0—90) ] s ' laminates as functions 

of the outer ply angle, 0. To investigate a and A^ in the region 

0 < 6 < 15 (and 75 < 0 < 90), the laminates [2/— 88] and [7.5/-82.5] were 

also considered. The ct values for these laminates are nearly equal to 

0.170, but the A^ values are much smaller than for the other laminates. This 

was expected because the stress is identically zero for [0/90] g and 

[90/0] laminates; thus, for 0 nearly equal to zero A should also be 
s 

very small. Dotted lines in figure 12 represent extrapolations to zero. 
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Figure 12 shows that both a and are symmetric about 0 — 45°. Also, 

the strength, A^, is maximum for the [15/— 75] g and [75/-l5] g laminates. 

[0/90] and [90/0] laminates.- For the [0/90] laminate, O z was used 
s s ® 

in the log-linear procedure because the was identically zero. As 

expected, the log-log plot in figure 13 shows a linear region with a negative 
slope. The power of the singularity was obtained from the least-square proce- 
dure as 0.202. This power is larger than the 0.17 calculated from other 
laminates but is only slightly beyond the range shown in figure 10. The com- 
puted power of 0.202 probably occurred because the a ^ singularity for the 
[0/90] laminate was not as dominant as the singularity analyzed for the 

other laminates. However, if the log-linear procedure were applied to O z 
stresses for smaller r values, the power may also approach 0.17. Such cal- 
culations for smaller r values would require further mesh refinement, which 
is beyond the scope of this paper. 

For the [90/0] laminate neither O nor 0 could be used in the 
log-linear procedure. The CJ ' was identically zero and the changed 

sign very near the free edge (see fig, 4(b)). Consequently, the singularity 
for the [90/0] laminate was not computed. However, as shown previously for 
other laminates, a laminate and its complementary laminate have the same a. 
Therefore, by induction, the singularities for [0/90] ^ and [90/0] | g should also 
be identical and a should be 0.17. 

In general, these results suggest that the power of the singularities for 

[0/ (6-90)] , 0 < 0 < 90, laminates with material properties used in this paper 

s 

is about 0.17. This power should be interpreted only as an approximate value 
because of the numerical accuracy of the log-linear procedure. 
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CONCLUDING REMARKS 

The stresses near a straight, free edge in composite laminates subjected 

to uniform axial strain were studied by a quasi-three-dimensional finite- 

element analysis. Seven different four-ply laminates were analyzed: [0/90] , 

s 

[ 15/— 7 5] , [30/-60] , £ ±45] , [60/-30] , [75/-15] , and [90/0] —belonging to 
s s s s s s 

the family [6/(6-90)] , 0 < 0 < 90. The laminates were idealized by eight- 

s 

noded isoparametric elements. Systematic convergence studies were made with 
the aid of three meshes to examine the interlaminar stresses near the free 
edge. 

Near the intersection of the interface and the free edge, the interlaminar 

normal stress 0 was tensile for [0/90]^ and [90/0] laminates and was com- 
z s s 

pressive for all other laminates considered. The interlaminar shear stress 

0 was much larger than the normal stress 0 for all laminates except for 

[0/90] and [90/0] laminates (for which a is identically zero) . Among 
s s xz 

the laminates considered, laminates [15/-75] and [75/— 15] developed the 

s s 

highest O 

XZ 

Convergence studies indicated that stress singularities probably exist at 
the intersection of the interface and the free edge, and the singularities do 
not exist at the midplane. A log-linear curve fitting procedure appears to 
confirm the existence of a singular point at the intersection of the interface 
and free edge. The shear stress data close to the singular point were used in 
the log-linear procedure to evaluate the power of the singularities. 

The present results suggest that the power of the singularity is the same for 
all orthotropic laminates in the family represented by [0/(0-9O)] g , 0 < 0 < 90 

for a given set of material properties. For these laminates, with material pro- 
perties of the graphite/ epoxy considered, the power of the singularity was about 


0.17. 
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